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ABSTRACT 



Aims. We present a study of the density profile of the remote M31 globular cluster B514, obtained from HST/ACS observations. 
Methods. Coupling the analysis of the distribution of the integrated light with star counts we are able to reliably follow the profile 
of the cluster out to r ~ 35", corresponding to ^ 130 pc. The profile is well fitted, out to ~ 15 core radii, by a King Model having 
C=1.65. With an estimated core radius = 0.38", this corresponds to a tidal radius of r, ~ 17" (~ 65 pc). The analysis of the light 
profile allows also the estimate of the ellipticity and position angle of the isophotes within r < 20". 

Results. We find that both the light and the star counts profiles show a departure from the best fit King model for r ^ 8" - as a 
surface brightness excess at large radii, and the star counts profile shows a clear break in correspondence of the estimated tidal radius. 
Both features are interpreted as the signature of the presence of extra tidal stars around the cluster. It is also shown that B514 has a 
half-light radius significantly larger than ordinary globular clusters of the same luminosity. In the My vs. logr/, plane, B514 lies in a 
region inhabited by peculiar clusters, like w Cen, Gl, NGC2419 and others, as well as by the nuclei of dwarf elliptical galaxies. 

Key words. Galaxies: individual: M 31 - Galaxies:star clusters - catalogs — Galaxies: Local Group 
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1. Introduction 

Until a couple of years ago we lack any knowledge of M3 1 clus- 
ters at large (projected) distances from the center of the galaxy; 
the farthest known cluster was Gl, at 7?^ ~ 35 kpc from the 
center of M31, while (a few) Galactic globulars are found out 
to galactocentric distance of Rqc ~ 120 kpc (10 at Rqc > 30 
kpc, according to Harris (1996)). Recent searches have identi- 
fied several new remote cluster in M3 1 (Huxor et al. (2004), 
Galleti et al. (2005), Bates et al. ( |2004 . Martin et al. (2006)). 
These studies seem to indicate that a relatively large number 
of globulars are still to be discovered in the extreme outskirts 
of M3 1 . A significant sample of distant globulars may provide 
extremely useful information on the early evolution of the halo 
of M3 1 . In particular, it has been suggested that bright globular 
clusters may be the remnants of disrupted nucleated dwarf galax- 
ies (Freeman & Bland-Hawthorn (2002), and references therein; 
see also Brodie & Strader (2006)): if this were the case the prob- 
ability of finding the observational fingerprints of these kind of 
phenomena is much higher at large distances from the center of 
parent galaxies, where substructures may survive for long times 
(Bullock & Johnston (2004)) and the overall stellar density is 
very low. Moreover, the structure and evolution of clusters orbit- 
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* Based on observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute (STScI), 
which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
** Table 5 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-t-A/ 

' In the present context we will call remote clusters those lying at 
galactocentric distances larger than 30 kpc. 



ing in very-low density environments is a very interesting topic 
in itself. 

In Galleti et al. ( (l2006bh . hereafter G06b) we have presented 
deep Hubble Space Telescope - Advanced Camera for Survey 
photometry of the recently discovered cluster B514 (Galleti et 
al. (2005), hereafter G05), lying at Rp ^ 55 kpc from the cen- 
ter of M3 1 . The derived Color Magnitude Diagram (CMD) re- 
vealed that B5 14 is a genuine old and metal poor globular cluster 
{{FefH} - -1.8, confirmed also by the spectroscopic estimate 
by G05). The cluster is very bright (My ^ -9.1) and appears 
quite extended, similarly to the brightest remote cluster of the 
Milky Way, i.e. NGC 2419. Here we present the analysis of the 
surface brightness distribution of B514 obtained from the same 
HST-ACS data. Coupling the surface brightness profile obtained 
from the integrated light - for the inner parts - to star counts 
in the outer region, and thanks to the extremely low level of 
background density in the field, we were able to identify an un- 
equivocal break in the outer profile of the clust er, ind icating the 
presence of extra- tidal st ars (see Johnston et al. (|1999[); C ombes. 
Leon & Meylan (1999); Grillmair et al. (119951: 1 19961) ; Leon, 
Meylan & Combes (2000)). Extra-tidal components and/or ex- 
tended tidal tails have be en obs erved in seve ral Ga lactic globu- 
lars (s ee Grillmair et al. (1 19951); L eon et al. (l2000h ; Testa et al. 
(120001) ; Oden kirchen et al. (l2003h ; Lee et al. (|2004 ; Bel okuroy 
et al. (120061) ; and references therein). Holland et al. (1 19971) 
and Barmby, Holland & Huchra (2002), hereafter BHH, found 
some M31 c lusters whose light profile exceeds the best fit King 
(Il962l;ll966h model in the outermost regions, and interpret ed this 
discrepancy as an extra-tidal component. Grillmair et al. (Il996l) 
found the same kind of discrepancy in three M31 clusters; they 
were able to follow the density profile of the clusters to signifi- 
cantly below the background level by coupling the light profile 
with the profile obtained from star-counts. By applying a simi- 
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lar technique, we are able to follow the profile of B514 out to 
r ~ 35". Moreover, we found that B514 has a half-light radius 
(r/,) larger (by ^ 15%) than typical globular clusters of the same 
luminosity, a characteristic shared by a few very peculiar sys- 
tems, like o ) Cen M54, Gl and NGC 2419 (see Mackey & va n 
den Bergh (|2005|) . hereafter MB05, and Hasegan et al. (|2005|) . 
for a thorough discussion). 

2. The surface brightness profile of B514 

Since the details of the observations and data reduction are re- 
ported in G06b, here we recall just a few essential elements, 
referring the interested reader to that paper. The cluster has 
been observed with the Wide Field Channel of the ACS. The 
WFC has a total field of view of 202" x 202" and a pixel scale 
of 0.049" pixel"'. The observational material is constituted by 
three F606W images (total f(,.vp = 2412 s) and three F814W im- 
ages (total tt,xp - 2419 s), and the associated combined (driz- 
zled) images. The cluster is placed in the center of one of the two 
ACS/WFC chips (Chip 2), while Chip 1 sample the field popu- 
lation (see G06b and below). If not otherwise stated, magnitudes 
are alw ays in the VEGAMAG scale as defined by Sirianni et al. 
(12005'). The reddening corrections are performed as in G06b, as- 
suming E{B - y) = 0.10; a distance modulus {m - M)o = 24.47 
is also assumed, after McConnachie et al. (2005), corresponding 
to Dq - 783 kpc. At this distance 1" correspond to 3.8 pc and 
one ACS/WFC pixel to 0. 19 pc. 

The extreme crowding conditions prevent the full resolu- 
tion into stars of the densest region of M31 globular clusters, 
even at the exquisite spatial resolution achieved by HST cam- 
eras. However, very accurate and well resolved surface bright- 
ness profiles can be obtained studying the distribution of thei r 
integrated light ( see, fo r example, Fusi Pecci et al. (|1994|) : 
Djorgovski et al. (l 2003l ); BHH). This technique, very success- 
ful in the bright inner regions of the clusters, is limited in the 
outermost part of the clusters, where the low luminosity density 
coming from cluster stars may be overwhelmed by the brightness 
of the background. In this regime star counts may be much more 
efficient, since cluster stars may be easily identified out to larg e 
radii, under favorable conditions (see Grillmair et al. (Il996l) ). 
The ACS field studied here lies more than three degrees apart 
from the center of M3 1 and it appears to have an exceedingly 
low density of background stars (see G06b and below). This al- 
lowed us to derive a reliable light profile, completely unaffected 
by incompleteness, out to r ~ 20", and to extend the analysis out 
to r ~ 35", in completely uncrowded regions, by counting stars 
having colors and magnitudes typical of the cluster population. 

2.1. The light profile 

The light profile has been obtained independently from the 
F606W and F814W drizzled images. A few heavily saturated 
foreground stars (all at r > 26" from the cluster center) have 
been excised from the images and replaced with the mean value 
of the surrounding background, to avoid contamination of the 
profile. 

The light profiles were derived using the XVISTA software, 
maintained by J. Holtzmaifl XVISTA iteratively resolves the 
profile by fitting ellipses to the observed light distribution un- 
til a stable solution is reached (see Fusi Pecci et al. (1994), 
Djorgovski et al. (l2003h . for examples of applications to M31 

^ See http : //astronomy . n msu . edu/holtz/xvista/index . html | 

for download and documentation 



clusters, and Lauer (Il985h . for a detailed description of the us- 
age). The code provides as output, the coordinates of the center, 
the surface brightness, the ellipticity (e = 1. - b/a), and the 
position angle (PA, in degrees, measured anti-clockwise from 
the North direction) of each fitted ellipse, as well as the total 
amount of light enclosed within each ellipse. The profiles are 
derived with a single pixel step. This resolution is appropriate 
for the innermost regions of the cluster (r < 2" - 3") where the 
light intensity is very high, while at larger radii provides a quite 
noisy profile. This problem is solved by getting an average of the 
whole profile over 10 px bins: this sampling ensure a satisfying 
level of noise out to r ~ 20". The background level is estimated 
as the average surface brightness in large (~ 100 x 100 px^) 
"empty" areas far away from the cluster The level of back- 
ground is very low (//feoew ~ 29.5 mag/arcsec^) and the final 
background-subtracted profiles were verified to be very robust 
to variations of the adopted background. The uncertainty on the 
position of the center is of order of ~ 2 px in the X and Y direc- 
tions: variations on the assumed position of this size doesn't af- 
fect significantly the derived profiles. The derived intensity pro- 
files can be directly converted into magnitude/arcsec"- units us- 
ing the VEGAMAG zero points of Sirianni et al. (2005). The 
derived F660W and F814W profiles are shown in Fig. [T] The 
"x"s in the innermost 2" are from the original 1 -pixel step pro- 
files, while the open circles with error bars are the averages in 
10 pixels bins. The profiles are quite smooth, well-behaved and 
very similar in shape, at least out to r ^ 8", i.e. the radius en- 
closing ^ 90% of the whole cluster light. Outside this radius 
both profiles sho w a ma rginal excess of light with respect to the 
best-fitting King (1962) model (see below), but the F814W pro- 
file appears more noisy, probably because of the larger weight 
associated to the contribution of individual bright RGB stars in 
this redder passband. The F606W light profile is reported in 
Tab. 3. The apparent integrated magnitudes in the VEGAMAG 
system mag,^f(,Q6w and mflg/ ^gHty were estimated by integrating 
the respective light profiles. 

We take the Half Width at Half Maximum (HWHM) of the 
profile as the core r adius (King (119621) . Spitzer (Il987h '). Once 
fixed this parameter we searched for the King's models provid- 
ing the best-fit to both profiles. In order to take in the proper 
account the ACS Point Spread Function (PSF) in measuring the 
light profile, the King models have been convolved with analytic 
F606W/F8 14W PSFs models, modeled on observed bright stars, 
as done by Barmby et al.(2007). All the comparisons between 
observed profiles and King's models presented in the follow- 
ing involve only PSF-convolved theoretical profiles. In Tab. 1 
we report both the observed and the de-convolved best-fit pa- 
rameters (see below). The former must be adopted when dealing 
with the observed profiles, while the latter must be used in the 
comparisons with other clusters. We note that the adoption of 
PSF-convolved profiles results in small changes 10%) in the 
cluster parameters, as usually occurs for extended M3 1 clusters 
like B514 (see Barmby et al.(2007)). 

The derived F660W and F8 14W profiles are shown in Fig.[Tl 
compared to PSF-convolved King's models with concentration 
parameter C - log(r,/rc) = 1.4, 1.5, 1.6, 1.7. If we limit to the 
most reliable region within r - 8" (encompassing ^ 19 core 
radii), the best fit is achieved with the C = 1.7 model for the 
F606W profile and C = 1 .6 for the F8 14W profile. 

We can gain some insight of the uncertainty associated with 
many observed/derived parameters by the comparison between 
the estimates obtained in the F606W and F814W profiles. For 
size parameters, (core radius r^, half-light radius /-/,, tidal radius 
r,, see King (il962,) ) we adopt the mean of the independent es- 
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Fig. 1. Surface brightness profiles (in mag/arcsec") of B514 in 
F606W (upper panel) and F814W (lower panel). The "x"s in the 
innermost region are from the 1 -pixel step profile, while the open 
circles are average values over 10 px bins. The curves are PSF- 
convolved King models of different concentration, from C=1.4 
to C=1.7, from left to right. The adopted core radius and central 
surface brightness are also reported. 



Fig. 2. Upper panel: color profile of B514; the continuous line 
is the overall mean color and the dotted lines encloses the +1 
standard deviation range. Middle panel: ellipticity profile. Lower 
panel: position angle profile. The latter two profiles are the aver- 
age of the profiles obtained from the F606W and F8 14W images. 
The meaning of the lines is the same as in the upper panel. 



timates obtained by the two profiles. For example, for the core 
radius we obtain rc - 0.38" from F606W and rc - 0.45" from 
F814W, and we adopt r,. = 0.42 + 0.03"; adopting C=1.65 we 
obtain a tidal radius r, ^ 18.8 + 2.5"; for the half-light radius 
we obtain r,, = 1.52" fromF606W and r/, = 1.73" fromF814W, 
and we adopt r/, = 1.6 + 0.2" (observed values, see Tab. 1). 

The half-light radii have been computed also by perform- 
ing aperture photometry on circular concentric annuli; this in- 
dependent procedure provided the same results obtained with 
XVISTA, indicating that the estimate of this parameter is very 
robust. Also the estimates of the apparent integrated magnitudes 
have been checked in this way, and the results obtained with dif- 
ferent methods are fully consistent. 

A summary of the measured structural parameters is pre- 
sented in Tab. 1, while Tab. 2 shows the derived parameters, i.e. 
those involving assumptions on distance and reddening and/or 
transformations to the standard Johnson-Kron-Cousins photo- 
metric system. The latter are achieved with the transformations 
presented in G06a. Note that if Sirianni et al.'s transformations 
are used instead, slightly brighter V magnitudes are obtained (by 
~ 0.06 mag), while the final I magnitudes are the same to within 
+0.02 mag. r, and surface brightness measures reported in 
Tab. 2 are derived from average de-convolved values. Hr,, is the 
mean surface brightness within r/,, while yu(0) is the central value 
of the surface brightness. 

In a very recent paper, Mackey et al (12007') reported on the 
CMDs and (observed) half-light radii and integrated magnitudes 
of eight M31 clusters, including B514 (their GC4), from inde- 
pendent ACS observations. The parameters of B5 14 obtained by 
this team are in go od agree ment with those presented by us here 
and in Galleti et al. (l2006bl) . 



2.1 .1 . Color profile, ellipticity and position angle 

In the upper panel of Fig. |2] we show the color profile of B5 14. 
The profile is remarkably constant in the inner 6", while it be- 
comes very noisy in the outer regions. It is clear that in the 
low-surface-brightness outer parts the contribution of individual 
RGB/BHB stars may be important to establish the local color. A 
general tendency toward redder colors in the outer regions can be 
noted, a phenomenon that has be en observed also in oth er clus - 
ters (see Djorgovs ki et a l. (II991I) . Djorgovski & Piotto ( Il993h . 
and Barmby et al. ( |2002| )). However, in the present case, we re- 
gard this trend as of marginal significance, given the uncertain- 
ties. 

The middle and lower panels of Fig.|2]display the ellipticity 
and position angle profiles, respectively. In both cases, once ver- 
ified that the profiles obtained from the F606W and F8 14W were 
fully consistent, we averaged the two. B514 has a mean elliptic- 
ity (e) = 0.19 ± 0.07, quite high for a globular cluster, but not 
extraordinary (MB05). A sizable enhancement of the ellipticity 
is apparent between r ~ 8" and r ~ 15". A twist of the isophotes 
seems to occur in the same radial range (lower panel, PA chang- 
ing by ~ 30°), suggesting a disturbed morphology in this range. 
The overall conclusion is that the cluster is rather elongated in 
the NW-SE direction, as is apparent by simply looking at the 
image (G06b). 

2.2. The profile from star counts 

For star counts we adopt the same catalog as G06b, including 
only well measured stars (see G06b, for details). Table 5 (online 
material) reports the photometry of the individual stars. In the 
upper panels of Fig. [3]we show the CMDs for the chip contain- 
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Table 1. Observed and de-convolved parameters 







— Observed — 




De-convolved 


Parameter 


F606W 


F814W 


Average 


Average 


[arcsec] 


0.38 


0.45 


0.42±0.03 


0.38 


rh [arcsec] 


1.52 


1.73 


1.6±0.2 


1.44 


r, [arcsec] 


19.0 


18.0 


18.8±2.5 


17.0 


C 


1.7 


1.6 


1.65 


1.65 


iu(0) [mag/arcsec-] 


16.41 


15.83 




16.33/15.74 " 


fir^ [mag/arcsec^] 


18.5 


17.9 




18.4/17.5 


mag, [V EG AM AG] 


15.48±0.06 


14.71±0.06 







" De-convolved F606W/F814W central surface brightnesses. 
De-convolved F606W/F814W half-light radius surface brightnesses. 



Table 2. Derived Parameters 



Parameter 


Average/ Adopted 


Mv 


-9.1+0.1 


M, 


-9.9 + 0.1 


l^viO) [mag/arcsec^] 


16.5 " 


/j/(0) [mag/arcsec^] 


15.7" 


yfv.r,, [mag/ai-csec^] 


18.4" 


fii^n, [mag/arcsec'] 


17.5 " 


r, [pc] 


1.4" 


n, [pc] 


5.4" 


(Ellipticity) 


0.19 ±0.07 


(Position Angle) [deg] 


28± 14 



" De-convolved quantities. 



ing the cluster (Chip 2) and for the chip presumably sampling the 
field population (Chip 1, see G06b). The reported contour is the 
filter we adopt to select likely cluster members on the CMD: it 
encloses Red Giant Branch (RGB) and Horizontal Branch (HB) 
stars having F814W < 25.5. The filter efficiently excludes ob- 
vious color outliers and faint stars whose membership can be 
uncertain. In the lower panels the X,Y map of the two samples 
- in their relative positions - is presented. The horizontal contin- 
uous line marks the boundary between the two chips. The stars 
selected by the filter are plotted as heavy dots. The larger cir- 
cle has a radius of 50" and is the largest circle that can be fully 
enclosed within one WFC chip. The following analysis is re- 
stricted only to filter-selected stars within this circle. The back- 
ground level of the stellar density is estimated from the whole 
50" circle in Chip 1 as pbkg = 0.0043 + 0.0013 stars/arcsec^ 
while we derive the cluster profile from selected stars in Chip 
2. It can't be excluded that cluster stars are present also in Chip 
1. However, estimating the background level in the X < 50", 
X > 150" regions of both chips (enclosed by dotted lines in the 
lower right panel of Fig. [3]), and in the Y < -50", Y > -50" 
regions of Chip 1 (separated by the long dashed horizontal line 
in the lower right panel of Fig. |3j, we found that the background 
is the same as that measured in the 50" circle, within < 2 - cr, 
ranging from 0.0024+0.0010 stars/arcsec^ to 0.0056+0.0015 
stars/arcsec^, and there is no discernible density gradient out- 
side r - 50" from the cluster center To have a more reliable 
estimate of the background level we would need observations 
of a larger (or more distant) field that, unfortunately, is lacking. 
However, this implies that, if anything, we are slightly overesti- 
mating the background. The possible associated bias would act 
against the detection of feeble extra-tidal components, hence it 



cannot be at the origin of the excess of surface brightness at large 
radii that is discussed below. 

The derived profile is shown in Fig. |4^. The profile is very 
extended: a level of 5-o- above the background is reached at r ^ 
31", while the (adopted) background level is reached at r ^ 35". 
This plot shows one of the main results of the present paper: 
not only the observed profile clearly extends much beyond the 
tidal radius derived from the light profiles but, above all, a clear 
change of slope is detected at r ~ 18", i.e. near r, itself. Note 
that the excess between r ^ 18" and r ~ 30" is many cr above 
the background, hence it is very significant, even if it encloses 
just a tiny fraction of the total cluster light. 

Given the extreme crowding conditions in the inner part of 
the cluster and the strong density gradient, it is expected that 
the completeness of the sample is subject to radial variations. 
Fig.|4]3 clearly illustrates the actual case by comparing the light 
profile (that is completely unaffected by incompleteness) and the 
background-subtracted star counts profile (reported in Tab. 4). 
A linear radial scale is adopted to provide a clearer comparison. 
This plot shows that incompleteness significantly affects star- 
counts for r ^ 6", becoming more and more important toward 
the center of the cluster, until, as said, it reaches 100% at r ^ 2". 
However an excellent match of the profiles can be achieved in 
the range 10" < r < 20" (i.e. to the end of the light profile). 
This clearly proofs that for r > 10" there is no more variation 
of the incompleteness with radius and, consequently, star counts 
provide a fair and fully reliable description of the real profile in 
the considered range. The vertical shift applied to the star counts 
profile to match the light profile automatically provides also the 
normalization constant to transform surface stellar densities into 
surface brightness. Therefore, the two profiles can be joined into 
one, covering the full < r < 35" range, as shown in Fig.|5] 

Is it possible that the excess component beyond r, and/or the 
observed change of slope can be due to sources unrelated to the 
cluster? This possibility is very hard to conceive, since (1) the 
adopted quality selections (G06a) and the CMD filter limit the 
analysis to relatively bright, well behaved stars that shouldn't 
suffer from any serious contamination, and (2) it is very hard 
to imagine a "field" population whose surface density decreases 
with distance from the cluster center. We must conclude that the 
detected surface density excess at large radius and the change 
of slope in the profile are genuine properties of the cluster. A 
change of slope in the outer regions of the surface density profile 
is generally interpreted as the signature of the presence of tidally 
stripped stars (see Combes et al. ( 1999), Johnston et al. ( 1999h : 
Leon et al. (2000), and references therein). For brevity, in the 
following we will refer to the stars beyond the break in the profile 
as to extra-tidal stars. 
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Fig. 3. Selection of likely cluster members on the CMD. Upper 
panels: CMDs of the stars in the ACS/WFC chip containing 
the cluster (left panel) and of those in the "empty" chip (right 
panel). The contour is the filter adopted to select star for fur- 
ther analysis; the dashed line marks the adopted magnitude limit 
(F814W< 25.5). Lower panels: maps of the selected (heavy 
dots) and unselected (but brighter than F814W= 25.5) stars of 
Chip 2 (above the continuous horizontal line) and Chip 1 (below 
the continuous horizontal line). Here we adopted a representa- 
tion that approximately displays the actual relative position of 
the two chips. The regions selected for the analysis are enclosed 
within the large circles (radius of 50"). This is the largest circle 
that can be enclosed in a single chip. The smaller circle in the 
lower-left panel has r - 20" and is indicative of the derived tidal 
limit of the cluster. The dotted vertical lines and the horizontal 
dashed line plotted in the lower right panel show the portions 
of the field that have been used to obtain estimates of the back- 
ground level, that were compared to that obtained in the 50" 
circle in Chip 1 (see text). 



3. Discussion 

A general prediction of incompleteness theoretical studies of 
tidal tails is that the slope of the surface brightness profile is 
different for bound stars and extra-tidal stars (see, for exam- 
ple. Combes et al. ( Il999l) . C99; Yim & Lee (2002); Johnston 
et al. (Il999h ; Montuori et al. (|2007|) ). Johnston et al. (i 19991) pre- 
dicts that the surface density of stars in the tails should decrease 
as r"', in agreement with most observations of ex tra-tidal stars 
around Galactic globulars (Grillmair et al. ('1995'); Leon et al. 
( 120001) ; Testa et al. (2000); Odenkirchen et al. (2003)). On the 
other hand C99 conclude that that such shallow slope is to be 
expected in the tidal tails at large distances from the parent clus- 
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Fig. 4. Panel (a): observed surface density profile obtained from 
star counts (grey squares). The continuous horizontal line marks 
the background level, the dotted lines marks the + 1 - cr levels 
around the background, and the dashed line marks the 5 - cr 
level above the background. The continuous vertical line marks 
the position of the tidal radius as derived from the fit of the light 
profile, assuming C= 1 .6, the long-dashed lines show the position 
of r, under the assumptions C=1.55 (left) and C=1.75 (right). 
Panel (b): Comparison and match between the F606W light pro- 
file (empty circles) and the (background-subtracted) star counts 
profile (grey squares). The scale of the horizontal axis has been 
kept linear to provide a clearer view of the superposition be- 
tween the two profiles. The dotted lines encloses the region in 
which the match between the two profile is excellent. 



ter - i.e. fully unbound independent tidal debris -, while in the 
vicinity of the cluster - i.e. immediately beyond the tidal radius - 
the density should decrease as r " with a = 3 or larger, and the 
involved stars cannot be considered as completely unbound. C99 
explain the discrepancy with the observed slopes as due to im- 
perfect subtraction of a very noisy background, typical in most 
Galactic cases. It is interesting to note that in the present case, 
where we deal with stars in the proximity of the cluster and the 
background is virtually non-existing, the density of extra-tidal 
component decreases significantly faster than r"', and it is com- 
patible with ps oc r^^ (the thick segment superposed to the outer 
profile in Fig.|5]l. 

In any case, the only possible alternative to explain the ob- 
served excess of stars and the change of slope in the profile at 
r, ^ 18" it is to postulate that the cluster is embedded in a 
very low surface brightness stellar system, that is, for instance 
an unknown dwarf galaxy. We consider this hypothesis as un- 
likely, since (a) the density of the detected extra-tidal compo- 
nent decreases with distance from the center of the cluster, that 
would be possible only if the cluster resides at the center of 
the hypothesized system,(b) the surface brightness of the extra- 
tidal component is t 26 mag/arcsec^, significantly lower than 
the typical central surface density of l ocal d warf spheroidals 
(Mateo a 998) ; but see Zucker et al. (|2004|) . for a counter- 
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Fig. 5. Total surface brightness profile obtained by joining the 
light (dark grey circles) and the star-counts (light grey squares) 
profiles. The light curve is the (PSF-convolved) King model best 
fitting the light profile. The dark segment in the outer part is a 

surface density oc r'^ power law. 



examp le in M31), and (c) the INT/WFC survey (Ibata et al. 
(1200 ll) ) reaches the surroundings of B514 (the cluster is near 
one edge of the survey field) and no particular substructure is 
detected in this region. The twisting and distortion of the outer 
isophotes noted in Sect. 2.1.1, above, also militates for the tidal 
origin of the observed structures. The detection of tidal tails on 
opposite sides with respect to the center of the cluster would 
have provided a conclusive argument in this sense (see Montuori 
et al. (120071) . and references therein). Unfortunately the extra- 
tidal over-density appear as statistically significant only when 
azimuth-averaged. We find just a small degree of anisotropy in 
the spatial distribution of extra-tidal stars: they appear to be lo- 
cated preferentially along the SE - NW direction, but the effect 
is not statistically significant. 

3.1. Bright clusters and tlie nuclei of dwarf ellipticals 

Extragalactic surveys performed with high spatial resolution 
cameras are revealing the existence of new kinds of s tellar sys- 
tems. The Extended Clusters (ECs) found by Huxor ( |2005|) are 
a typical example, but other kinds of extended cluste rs hav e 
been identified in more distant galaxies (see Peng et al. (2006|), 
for a thorough review and discussion). HiUcer et al. (1999) and 
Drinkwater et al. ( l2000h recently discovered a new kind of dwarf 
galaxies inhabiting galaxy clusters that are slightly more lumi- 
nous than the brightest globular clusters and significantly more 
compact than any dwarf galaxy, the Ultra Compact Dwarf galax- 
ies (UCD). It has been suggested that such systems are the dense 
remnant of tidally harassed nucleated galaxies (Drinkwater at al. 
( l2003h ). In this context a new perspective on the possible rela- 
tions and differences among these various system is gradually 
emerging, and the comparison among the structural properties 
may provide interesting insights in this sense. 



Table 3. Background-subtractedF606W Surface Brightness pro- 
file from the integrated light. 



r 






arcsec 


mag/arcsec^ 


mag/arcsec^ 


0.24 


16.7 


0.1 


0.73 


17.9 


0.1 


1.23 


18.8 


0.2 


1.71 


19.7 


0.1 


2.20 


20.2 


0.2 


2.69 


20.5 


0.2 


3.18 


21.2 


0.2 


3.67 


21.6 


0.2 


4.16 


22.1 


0.2 


4.65 


22.2 


0.2 


5.14 


22.5 


0.4 


5.63 


22.8 


0.2 


6.12 


23.1 


0.3 


6.61 


23.5 


0.3 


7.11 


23.1 


0.6 


7.59 


23.5 


0.4 


8.09 


24.0 


0.2 


8.57 


24.3 


0.4 


9.06 


23.9 


0.3 


9.55 


23.9 


0.4 


10.05 


24.1 


0.4 


10.53 


24.8 


0.5 


11.02 


24.6 


0.3 


11.51 


24.8 


0.3 


12.01 


24.5 


0.8 


12.49 


25.0 


0.5 


12.98 


25.4 


0.4 


13.47 


25.8 


0.4 


13.96 


26.0 


0.3 


14.45 


25.7 


0.4 


14.94 


26.2 


0.6 


15.43 


26.1 


0.6 


15.92 


26.2 


0.6 


16.41 


26.2 


0.6 


16.90 


26.7 


0.6 


17.40 


26.7 


0.7 


17.89 


26.7 


0.6 


18.38 


26.6 


1.4 


18.86 


27.0 


0.5 


19.35 


27.0 


0.6 



In this line, MB05 compared globular clusters from different 
galaxi es in th e My vs. logr/, diagnostic plane (see also Hasegan 
et al. (l2005l) ). They found that the bulk of globulars lies be- 
low the line logr;,=0.25Mv+2.95. All the objects that are found 
above this threshold, namely to Cen, NGC 2419, M54 and Gl, 
are very bright and anomalous clusters: all of them were pre- 
viously indicated as possible remnants of disrupted nucleated 
dwarf galaxies (Freeman & Bland-Hawthorn ( |2002|) . The My vs. 
logri, plane has a distinct advantage with respect to other simi- 
lar dia gnostic planes (as, for instance. My vs yUy(O), Kormendy 
('1985')), since r/, is a quite easy-to-measure and reddening inde- 
pendent quantity. 

In Fig.|6]we show the position of B514 and other interesting 
systems in the My vs. logr/, plane. Filled circles are Galactic 
globular clusters, from MB05; grey crosses are M31 globu- 
lars, our own estimates from HST imaged; open circles are 

^ These estimates have been obtained with a profile analysis strictly 
homogeneous to that performed here for B514; the results have been 
compared with the independent estimates by BHH: the agreement be- 
tween the observed parameters from the two sources is very good. 
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Table 4. Background-subtractedF606W Surface Brightness pro- 
file from star counts. 



r 

arc sec 


mag/arcsec 


mag/arcsec^ 


10.41 


24.7 


0.1 


11.64 


25.0 


0.1 


12.86 


25.3 


0.1 


14.09 


25.8 


0.1 


15.31 


26.2 


0.1 


16.54 


26.6 


0.2 


17.76 


27.0 


0.2 


18.99 


27.1 


0.2 


20.21 


27.2 


0.2 


21.44 


27.4 


0.2 


22.66 


27.5 


0.2 


23.89 


27.8 


0.2 


25.11 


28.0 


0.2 


26.34 


28.1 


0.2 


27.56 


28.4 


0.2 


28.79 


28.6 


0.2 


30.01 


29.1 


0.3 


31.24 


29.5 


0.3 


32.46 


29.6 


0.3 


33.69 


29.9 


0.3 


34.91 


29.8 


0.3 



" The star counts profile has been converted into surface bright- 
ness units (mag/arcsec-) with the relation yUf^sosiv = -2.5log(ps)+ 
23.20, according to the normalization shown in Fig.lH 

the M31 clusters of Mackey et al. (l2007b . seven of them hav- 
ing Rp > 30 kpc (in the following we will refer to these clus- 
ters as M-GCl, M-GC2, M-GCIO); asterisks are ECs from 
Mackey et al. (,2006); open stars a re UCDs in the Fornax clus- 
ter from De Propris et al. (l2005h . The continuous line is the 
Mackey & van den Bergh threshold mentioned above. Apart of 
B514, the real novelty of Fig. ^ with respect to previo us ver- 
sions (MB05; Huxor et al. (120051) : Belokurov et al. (l2007l) ) is that 
for the first time it is possible to report also the position of the 
nuclei of dwarf nucleated ellipticals (red triangles). Cote et al. 
(12006*) measured half-light radii and gAB, ZaelI integrated mag- 
nitudes of the nuclei of several dwarf ellipticals in Virgo from 
deep ACS/WFC images. We converted g^B magnitudes into V 
magnitudes with the transformation: 

V - gAB - Q3l{gAB - Zab) (r.m.s. - 0.1 t7iag) 

that we have obtained from 166 bright stars of the cluster 
NGC2419 that are in common between the gAB, ^ab photome- 
try we obtained from archive ACS/WFC data and the secondary 
standards provided by Stetson (2005) for this cluster Then, we 
converted the integrated V magnitudes into absolute magnitudes 
by adoptin g the reddening and distance modulus provided by 
Cote et al. [20061) . 

There are a number of very interesting indications emerging 
from Fig.|6] 

1. In addition to the clusters already noted by MB05, i.e. M54, 
CO Cen, NGC 2419 and Gl, there are a few other bright M31 
globulars lying above the " ordinary globular cluster" thresh- 
old: B514, M-GC3, M-GC5, G76, G280, G219, and G302. 
All of these systems are classified as globular clusters but 
have half-light radii significantly larger than those of typical 
genuine globular clusters of the same metallicity. 

" gAB and Zab are the ACS/WFC F475W, F850LP pas sbands , respec- 
tively, calibrated in the ABMAG system; see Cote et al. ll2006h . 




Fig. 6. Absolute integrated V magnitudes vs. logarithm of the 
half light radius in pc for various low-luminosity stellar systems. 
The plotted line is the threshold for ordinary globular clusters in 
this plane as defined by MB05, log r/, = 0.25/Wi/-i-2.95. Filled 
circles are Galactic globular clusters, from MB05; grey crosses 
are M31 globulars, our own estimates from HST images; aster- 
isks are ECs from Mackey et al^2006); open circles are outer 
M3 1 GCs from Mackey et al. (i2007b ; open stars are UCDs in the 
Fornax cluster from De Propris et al. ( l2005i) ; red triangles are nu- 
clei ofdwarf elliptical galaxies in the Virgo cluster from Cote et 
al. ('2006'); the empty circles are the M31 clusters of Mackey et 
al. (2007);the grey square is B514. The globular clusters lying 
above the line are labeled with their names. 



2. The distribution of nuclei nicely overlaps the position of 
these anomalous clusters. In particular, the nuclei of dwarf 
ellipticals appear to join the brightest globulars to the anoma- 
lous "above-threshold" clusters and to the UCD galaxies, 
which have been also interpreted as th e nucl ear remnants 
of shredded galaxies (Drinkwater et al. (|2003|) ). This is the 
first time that a clear connection between the structure of 
nuclei, UCDs and anomalous clusters is directly established 
by comparing sizable samples. It would be of great interest 
to extend the comparison to velocity dispersions, but unfor- 
tunately these quantities are not available for th e sam ple of 
nuclei considered here (but see Hasegan et al. (l2005b ) and 
for most of M3 1 globulars, including B5 14 . 

3. ECs, on the other hand, seem to have a diff'erent nature, 
somehow intermediate between globular and open star clus- 
ters (Peng et al. (l2006h ). 

4. The remote clusters B514 {Rp ^ 55 kpc) and M-GC5 (Rp - 
78 kpc) lie above the logr/, - My threshold, while M-GCl 
(Rp ^ 46 kpc) and M-GCl (Rp - 100 kpc) are located well 
below the threshold, fully immersed within the distribution 
of ordinary globular clusters. This demonstrates that a large 
half-light radius is not a distinctive characteristic of remote 
clusters. 

5. It is interesting to note that the faint Galactic satellites re- 
cently discovered by various SDSS teams (see Belokurov et 
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al. (|2007), and references therein) lie outside the limits of 
Fig. |6] all of them having logr/, > 1.6 and My > -8, just 
above ECs. Ordinary dwarf sp heroid als have similar sizes 
but they have My < -8 (Mateo (Il998l) ). 

It is interesting to have a closer look to the anomalous 
"above-threshold" clusters. a> Cen and Gl clearly host stars of 
different chemical composition (and, presumably, age), hence - 
unlike classical globulars - t hey we re able to sustain c hemica l 
evolution (see, SoUima et al. ( l2005h : Bekki & Freeman ( |2003|) : 
Meylan et al. ('2001'), and references therein). Both clusters are 
quite elliptical in shape, similar to B514. M54 resides wit hin the 
nucleus of the Sagittarius dwarf spheroidal (Monaco et al. (l2005h 
and there are indications of a small metallicity spread among its 
stars (Sarajedini & Layden (Il995h ). G7 6 is a relatively metal 
poor cluster ([Fe/H]^ -1.3, Rich et al. (200$), hereafter R05) 
that is projected onto a very d ense star forming region in the 
disc of M31 (Bellazzini et al. (l2003l) ). The extreme crowding 
conditions in this region prevents a clear interpretation of the 
wide RGB shown in its CMD as obtained by R05, but this fea- 
ture clearly leaves room for a possible metallicity spread. G280 
is a quite metal rich ([Fe/H]^ -0.5, R05) and bright cluster; like 
G76 it is projected onto a high density background. G302 is a 
metal poor cluster ([Fe/H]^ -1.7, R05) with a blue horizontal 
branch. Its CMD (R05) is quite clean and it suggest that G302 is 
a normal single -population cluster G219 is a metal poor cluster 
( [Fe/H] ^-1.9, R05) lo cated at ~ 20 kpc from the center of M3 1 . 
Bellazzini et al. (I2003h noted that G21 9 is pro jected onto the gi- 
ant stream discovered by Ibata et al. (1200 Ih . It is remarkable, 
in the present context, that both G302 and G219 are among the 
rare M31 clusters with a detection of extra-tidal stars, as B514 
(Holland et al. (1997); GriUmair et al. (1996)). NGC2419 is a re- 
mote metal-poor Galactic gl obular whose CMD is very similar 
to that of B514 (Harris et al. (1991)). As B514 it show no obvi- 
ous sign of a metallicity spread, but this may be very difficult to 
find out from the CMD of very metal-poor clusters (see MB05). 
It has to be noted that NGC2419 has a half-light radius as large 
as those of the largest dE nuclei and UCD galaxies. 

In summary, B514 fully lies into the region of the My vs. 
logr/, plane that appear forbidden to ordinary clusters, in com- 
pany of a few other clusters, many of which present some kind 
of peculiarity. The fact that the same region hosts also the nu- 
clei of dwarf ellipticals provide support to the hypothesis that 
these bright extended clusters can be th e rem nant of disrupted 
galaxies (Freeman & Bland-Hawthorne (|2002|) ; MB05, and ref- 
erences therein; see Brodie & Huchra for a thorough discus- 
sion on the role of globular clusters within a cosmological con- 
text). On the other hand, it is possible that ordinary clusters 
are allowed to attain extended structures if their orbit never 
drives them in the inner part of the parent galaxy, as it may be 
the case of NGC 2419 and B514 (see van den Bergh, Morbey 
& Pazder (1199 lb ). However, the coexistence of ordinary and 
"above-threshold" clusters in the outermost regions of M3 1 (see 
point 4, above) does not support this hypothesis. 

In any case the present study confirms that globular clus- 
ters in the outskirts of M31 may reveal many interesting fea- 
tures, thanks also to the favorable observing conditions (as for 
instance, the very low background). Remote clusters are rare in 
the Milky Way, and they are typically quite faint, with the only 
exception of NGC2419 (see Mackey et al. (2007) and Galleti et 
al. (2007)). Hence, M31 may provide the opportunity of a sys- 
tematic study of a kind of stellar system that is very rare in the 
Galaxy. 



Using the same selection criteria that lead us to the discov- 
ery of B514 (G05) we have selected a conspicuous number of 
candidate clusters at large distances from the center of M3 1 ; at 
present we have obtained spectroscopic confirmation (as in G06) 
for fo ur of t hem, all being nearly as luminous as B514 (Galleti 
et al. (|2007|) ). The follow-up of these remote clusters, as well as 
of those discovered by other teams, may hopefully open a new 
window in the study of the M3 1 system and of bright clusters as 
a whole. 
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